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ARTICLE  INFO  ABSTRACT 


The  Sea-viewing  Wide  Field-of-View  Sensor  (SeaWiFS)  derived  diffuse  light  attenuation  along  the  Louisiana 
continental  shelf  (LCS)  was  examined  at  monthly  scales  from  1998  to  2007  to  characterize  temporal  and  spa¬ 
tial  patterns,  and  responsible  physical  forcing  conditions.  The  SeaWiFS  diffuse  light  attenuation  ranged  from 
0.10  to  2.64  m  \  Stepwise  multiple  linear  regression  analysis  suggested  that  spatial  and  temporal  patterns 
in  diffuse  light  attenuation  were  influenced  by  wind  speed,  nutrient  loading,  and  river  discharge  from  the 
Mississippi  and  Atchafalaya  River  Basin.  SeaWiFS  daily  integrated  surface  photosynthetically  active  radiation 
(PAR,  400-700  nm)  and  diffuse  light  attenuation  were  used  to  calculate  the  absolute  PAR  and  percentage  of 
surface  PAR  that  reached  the  sediment  water  interface  (SWI)  on  the  LCS.  Large  portions  of  the  LCS  were  eu- 
photic  to  the  SWI  especially  during  April  and  May.  This  finding  implied  that  significant  primary  production 
was  possible  beneath  the  pycnocline  during  spring  and  early  summer.  In  addition,  this  study  was  the  first 
to  demonstrate  that  the  euphotic  depth  was  correlated  to  the  depth  at  which  the  water  column  turned  hyp¬ 
oxic  on  the  LCS.  The  development  of  hypoxic  waters  may  be  influenced  by  decreased  light  availability  below 
the  pycnocline  in  addition  to  aforementioned  physical  forcing 
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1.  Introduction 

The  Louisiana  continental  shelf  (LCS)  is  strongly  influenced  by 
the  discharge  of  the  Mississippi  and  Atchafalaya  River  Basin  (MARB). 
Nutrient-rich  water  is  delivered  to  the  LCS  through  the  MARB,  with 
maximum  concentrations  occurring  in  the  spring.  Two  thirds  of  the 
Mississippi  River  flow  enters  the  Gulf  of  Mexico  through  the  Mississippi 
River  bird-foot  delta,  whereas  the  remaining  third  is  delivered  from  the 
Atchafalaya  River  (Dale  et  al.,  2010).  Historical  records  indicate  that 
nutrient  concentrations  in  the  Mississippi  River  increased  approximate¬ 
ly  2-fold  during  a  25-year  period  beginning  in  the  mid-1950s  (Goolsby 
et  al.,  200 1 ).  However,  there  was  no  significant  change  in  river  nutrient 
concentrations  from  1979  to  2007  (Greene  et  at,  2009). 

River  nutrient  loading  increases  primary  productivity,  phytoplank¬ 
ton  biomass  and  therefore  organic  matter  within  the  river  plume  and 
adjacent  coastal  waters  (Lohrenz  et  al.,  1997,  2008;  Turner  &  Rabalais, 
1991).  The  effects  of  nutrient  enhanced  productivity  across  the  shelf  re¬ 
gion  have  been  inferred  from  modeling  (Bierman  et  al.,  1994;  Justic  et 
al.,  1996;  Scavia  et  al.,  2003)  and  by  sediment  geochemical  measures 
(Eadie  et  al.,  1994,  Rabalais  et  al.,  2004,  2007;  Turner  et  al.,  2004).  The 
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transport  and  downward  flux  of  organic  matter  to  the  sub- pycnocline 
waters  is  believed  to  fuel  microbial  metabolism  that  consumes  oxygen, 
thus  contributing  to  hypoxic  conditions  below  the  pycnocline  (Redalje 
et  al.,  1 994;  Wysocki  et  al.,  2006). This  net  oxygen  consumption  contrib¬ 
utes  to  the  largest  zone  of  oxygen-depletion  in  the  United  States  coastal 
waters,  averaging  1 3,800  km1 2  between  1 985  and  2008  based  on  annual 
July  surveys  (http://gulfhypoxia.net). 

The  net  consumption  of  oxygen  in  below-pycnocline  waters  and  the 
sediments  is  the  primary  cause  of  hypoxia.  Potential  mechanisms  off¬ 
setting  oxygen  consumption  below  the  pycnocline  are  processes  such 
as  lateral  and  vertical  oxygen  flux,  and  in  situ  oxygen  production  via 
light  driven  photosynthesis  (Justic  et  al,  1996;  Lehrter  et  al.,  2009). 
Light  below  the  pycnocline  can  support  both  phytoplankton  and  the 
microphytobenthic  ( MPB)  net  production.  MPB  assemblages  can  persist 
at  light  levels  as  low  as  0.1  %  surface  irradiance  (Cahoon,  1 999),  whereas 
phytoplankton  typically  become  light-limited  at  1%  of  surface  irradi¬ 
ance  (Kirk,  1 994).  While  diffuse  light  attenuation  is  high  within  the  Mis¬ 
sissippi  River  plume  (Lohrenz  et  al.,  1999),  it  decreases  with  increasing 
salinity  such  that  large  portions  of  the  LCS,  where  hypoxia  occurs,  are 
euphotic  throughout  the  water  column  and  to  the  sediment  water  in¬ 
terface  (SWI)  (Lehrter  et  al.,  2009).  The  Bierman  et  al.’s  ( 1994)  deter¬ 
ministic  mass  balance  model  suggests  that  light  penetrating  below  the 
pycnocline  is  an  important  factor  influencing  organic  matter  production 
and  dissolved  oxygen  concentrations.  Dortch  et  al.  (1994)  estimated 


BA  Schaeffer  et  at  /  Remote  Sensing  of  Environment  115  (2011)  3748-3757 


3749 


that  photosynthesis  at  the  sediment- water  interface  contributed  up  to 
57%  of  the  lower  water  column  photosynthesis.  In  this  respect,  an  in¬ 
creased  understanding  of  light  availability  on  the  LCS  is  critical  to  un¬ 
derstanding  the  relative  contributions  of  water  column  and  SWI 
photosynthesis  to  the  supply  of  carbon  and  oxygen  in  the  hypoxic  zone. 

In  this  paper  we  examine  shelf-wide  synoptic  estimates  of  diffuse 
light  attenuation,  absolute  daily  integrated  photosynthetically  active 
radiation  (PAR,  400-700  nm),  and  the  percentage  of  daily  integrated 
PAR  at  the  SWI  using  satellite  observations.  We  characterize  the  tim¬ 
ing,  location,  and  physical  forcing  impacts  on  diffuse  light  attenuation 
and  PAR  reaching  the  SWI.  While  surface  PAR  and  diffuse  light  atten¬ 
uation  on  the  LCS  are  routinely  measured  through  a  variety  of  ship- 
based  surveys,  these  observations  are  typically  limited  to  a  few  sta¬ 
tions  per  day  spatially  and  only  a  few  weeks  at  a  time  temporally.  In 
addition,  field  surveys  in  the  northern  Gulf  of  Mexico  are  biased  to¬ 
wards  the  spring  and  summer  months  when  hypoxia  is  prevalent 
Therefore  field  observations  are  inadequate  to  fully  describe  water 
column  diffuse  light  attenuation  across  the  entire  LCS  over  the  entire 
annual  cycle.  In  contrast,  satellite-derived  data  provide  synoptic  cov¬ 
erage,  once  properly  validated  against  shipboard  or  buoy  observa 
tions.  The  goals  of  this  study  are:  (1)  validate  the  Sea-viewing  Wide 
Field-of-View  Sensor  (SeaWiFS)  daily  integrated  surface  PAR,  (2)  ad¬ 
just  the  SeaWiFS  Kd490  diffuse  light  attenuation  coefficient  product  to 
represent  KdPAR  with  field  observations  of  KdPAR,  (3)  estimate  the 
spatial  and  temporal  patterns  of  absolute  and  percentage  surface 
PAR  at  the  SWI,  (4)  interpret  how  the  euphoric  depth  relates  to  hyp¬ 
oxia,  and  (5)  examine  relationships  between  SeaWiFS  diffuse  light  at¬ 
tenuation  and  wind  speed,  MARB  discharge,  and  nutrient  loading 
from  1998  to  2007. 

2.  Methods 

Field  observations  were  used  from  EPA  Gulf  Ecology  Division 
shelf-wide  surveys  (Fig.  1 ),  Louisiana  Universities  Marine  Consortium 


(LUMCON)  July  surveys  (Fig.  1,  unpublished  data,  NOAA  NODC,  www. 
nodc.noaa.gov),  and  the  LUMCON  Lake  Pontchartrain  station  (Fig.  1). 

2.1.  EPA  Gulf  Ecology  Division  field  observations 

Seven  research  cruises  were  conducted  from  2005  to  2007  on  the 
LCS  where  hypoxia  routinely  occurs  (Fig.  1 ).  Three  cruises  were  con¬ 
ducted  in  the  spring  (March  2005,  April  2006,  and  April  2007)  and  4 
cruises  were  conducted  during  summer  months  (September  2005, 
June  2006,  September  2006,  and  August  2007).  The  March  2005 
cruise  was  aboard  the  R/V  Longhorn  and  the  remaining  cruises  were 
aboard  EPA’s  Ocean  Survey  Vessel  Bold.  With  the  exception  of  the 
September  2005  cruise,  a  Sea-Bird  911  CTD  (Sea-Bird  Electronics, 
Bellvue,  WA,  USA)  was  deployed  to  measure  vertical  profiles  of  tem¬ 
perature,  salinity,  dissolved  oxygen,  and  photosynthetically  active  ra¬ 
diation  (PAR,  Biospherical  Instruments,  QSP-200,  San  Diego,  CA,  USA). 
A  Sea-Bird  43  dissolved  oxygen  sensor  on  the  CTD  was  used  to  esti¬ 
mate  the  depth  at  which  the  water  column  turned  hypoxic  (defined 
as  <2  mg  02  L“  \  1 .4  ml  02  L”  \  or  -62  02  pmol  kg  1 ).  Secchi  depth 
was  recorded  at  each  station  occupied  during  daylight  hours. 

Surface  PAR  (400-700  nm)  (Eo)  was  measured  continuously  with 
a  2tt  U-COR  U-190  irradiance  sensor  and  averaged  at  15  min  inter¬ 
vals.  The  diffuse  light  attenuation  coefficient  (KdPAR,  m  ’)  deter¬ 
mined  from  the  CTD  was  calculated  as  the  slope  of  In(E2/Eo)  versus 
depth  (z)  where  E*  was  the  PAR  at  depth  z  and  Eo  was  the  PAR  mea¬ 
sured  at  the  surface.  The  reciprocal  of  Secchi  depth  (SD)  was 
regressed  against  KdPAR  at  Gulf  Ecology  Division  sites.  The  regression 
coefficients  were  then  used  to  estimate  KdPAR  at  all  LUMCON  sites. 

22.  LUMCON  field  observations 

Secchi  depth  was  used  to  estimate  diffuse  light  attenuation  coeffi¬ 
cients  during  LUMCON  surveys  because  PAR  measurements  were  not 
reported  at  LUMCON  sites.  SeaWiFS  surface  daily  integrated  PAR 
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Fig.  1.  Map  of  study  region  in  the  northern  Culf  of  Mexico.  Locations  of  all  EPA  Gulf  Ecology  Division  (GED)  sample  stations  are  indicated  in  open  circles.  LUMCON  sample  stations 
are  indicated  with  asterisks,  and  the  LUMCON  Lake  Pontchartrain  monitoring  station  is  indicated  as  a  diamond. 
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(Frouin  et  al.(  2003)  was  validated  against  the  LUMCON  Lake  Pontch- 
artrain  environmental  monitoring  station  (Fig.  1)  retrieved  on  Janu¬ 
ary  2009  from  the  archive  data  site  at  http://weather.lumcon.edu/ 
archivedata/select  asp.  The  monitoring  station  PAR  sensor  was  cali¬ 
brated  every  2  years  and  serviced  every  other  month  with  periods 
of  monthly  servicing.  Instantaneous  PAR  was  averaged  hourly  at  the 
station  and  integrated  throughout  each  day  from  2001  to  2007 
using  Interactive  Data  Language  (IDL,  ITT  VIS).  Data  for  daily  integra¬ 
tion  from  2004  were  not  available  at  the  Lake  Pontchartrain  station. 
LUMCON  dissolved  oxygen  measurements  were  made  with  a  Hydro- 
lab  Surveyor  3  (unpublished  data,  NOAA  NODC)  and  was  used  to  es¬ 
timate  the  depth  at  which  the  water  column  turned  hypoxic. 

23.  SeaWiFS  sate//ite  observations 

Satellite  ocean  color  data  were  obtained  from  the  National  Aero¬ 
nautics  and  Space  Administration  (NASA)  Ocean  Color  Web  (Feldman 
&  McClain,  2008).  SeaWiFS  provided  daily  images  with  pixels  having 
a  nominal  1.1  km  spatial  resolution  when  not  obstructed  by  cloud 
cover.  Atmospheric  correction  algorithms  applied  by  NASA’s  fifth 
(5.2)  SeaWiFS  reprocessing  included  the  standard  near-infrared 
(NIR)  correction  for  coastal  waters  (Arnone  et  al.,  1998;  Stumpf  et 
al.,  2003).  SeaWiFS  Data  Analysis  System  (SeaDAS)  version  5.3 
(Baith  et  al..  2001)  was  used  to  process  data  from  SeaWiFS  level-1 
to  level-2  one-day  images  for  comparison  against  field  observation 
data.  A  3  x  3  pixel  box  was  used  for  each  corresponding  station  lati¬ 
tude  and  longitude  coordinate  as  described  in  Bailey  and  Werdell 
(2006).  Satellite  match-ups  for  validation  of  PAR  and  diffuse  light  at¬ 
tenuation  were  evaluated  with  a  Model  II  geometric  mean  linear  re¬ 
gression  (Laws  &  Archie,  1981)  using  only  in  situ  data  collected 
between  10:00  and  14:00  local  time  on  the  same  day  as  the  satellite 
overpass.  Monthly  composites  were  temporal  averages  of  flag-free 
pixels  in  each  month  from  January  1998  through  December  2007. 
Monthly  composites  were  processed  with  the  Naval  Research  Labora¬ 
tory’s  Automated  Processing  System  (APS)  (Martinolich,  2005).  Sea¬ 
WiFS  imagery  covered  between  27.15  to  30.40°  N  and  88.22  to 
94.45°  W. 

Kd49o  (Lee  et  al.,  2005)  was  adjusted  to  represent  KdPAR  using  the 
slope  and  intercept  of  an  empirical  regression  between  KdPAR  field 
observations  and  the  SeaWiFS  Kd4qo  product.  The  exponential  light 
extinction  equation 

Ez  -  Eo*e,-Kd*z) 

was  used  to  calculate  the  absolute  and  percentage  of  surface  PAR  at 
the  SWI,  and  the  depth  of  the  euphotic  layer,  where  Ez  was  the  calcu¬ 
lated  light  at  the  SWI,  Eo  was  SeaWiFS  daily  integrated  surface  PAR, 
Kd  was  the  SeaWiFS  Kd4c*>  PAR,  and  z  was  the  depth  of  the  SWI. 

2.4.  QuickScat  wind 

Wind  speed  data  were  obtained  from  the  Physical  Oceanography 
Distribution  Active  Archive  Center  (PO.DAAC)  Ocean  Earth  Science  In¬ 
formation  Partner  (ESIP)  Tool  (POET)  with  the  Quick  Scatterometer 
(QuickScat)  from  27.15  to  30.40’  N  and  88.22  to  94.45°  W.  QuickScat 
measures  wind  speed  of  3  to  20  m  s“ 1  with  an  estimated  accuracy  of 
2  m  s-1  and  vector  resolution  of  25  km  (JPL,  2001).  Monthly  compos¬ 
ites  were  compiled  from  January  2000  through  December  2007. 

2.5.  Bathymetry,  river  discharge,  and  nutrient  loads 

Bathymetry  data  were  obtained  from  the  National  Geophysical 
Data  Center  (NGDC  Coastal  Relief  model,  retrieved  on  May  2009 
from  www.ngdc.noaa.gov/mgg/coastaI/coastal.html).  The  native  res¬ 
olution  was  3  arc  sec  (-0.1  km),  so  the  grid  was  resized  to  match 
1  km  resolution  of  the  satellite  imagery.  The  offshore  boundary  for 


the  domain  was  the  70  m  isobath.  Monthly  average  nver  discharge 
and  nutrient  loads  were  obtained  from  the  USGS  (http://toxics.usgs. 
gov/hypoxia/mississippi/nutrienc  flux  yield  est.html). 

A  stepwise  multiple  linear  regression  analysis  was  used  to  deter¬ 
mine  if  river  discharge,  nutrient  concentration,  and  wind  speed 
were  important  parameters  in  Kd  interannual  variability.  The  analysis 
equation  was: 

Kd  =  Wind  f  Q  +  Q.l  +  Q.2  +  NOx  f  NOx.l  4  NOx.2  4  Q:NOx 
4Q.l:N0x.l  4Q.2:NOx.2 

where  Q  was  the  river  discharge,  Q.l  and  Q2  were  the  one  and  two- 
month  lagged  discharges,  NOx  was  the  nitrate  concentration  in  the 
river,  NOx.1  and  NOx.  2  were  the  one  and  two- month  lagged  concen¬ 
trations,  and  Q:NOx  terms  were  the  interaction  terms. 

Four  alternate  future  scenarios  were  used  to  predict  change  in  dif¬ 
fuse  attenuation  coefficients  within  each  region  from  the  stepwise 
multiple  linear  regression  analysis.  Scenarios  were  selected  from 
Greene  et  al.  (2009)  to  forecast  changes  based  on  suggested  nutrient 
management  strategies  (Mississippi  River/Gulf  of  Mexico  Watershed 
Nutrient  Task  Force,  2001, 2008)  and  global  climate  change  scenarios 
(Miller  &  Russell,  1992;  Nijssen  et  al.,  2001 ;  Wolock  &  McCabe,  1999). 
Scenarios  included  (1)  30%  and  (2)  50%  reduction  in  May  NOx  con¬ 
centration  with  no  change  in  May  discharge,  (3)  20%  increase  in 
May  NOx,  and  (4)  20%  increase  in  May  discharge. 

3.  Results 

3.1.  Algorithms 

There  were  1257  paired  observations  from  the  LUMCON  Lake 
Pontchartrain  environmental  monitoring  station  and  the  SeaWiFS  im¬ 
agery  taken  from  2001  through  2007.  The  paired  observations  were 
randomly  divided  so  half  of  the  data  were  used  to  develop  regression 
coefficients  and  the  other  half  of  the  data  were  used  to  validate  the 
regression  analysis.  There  was  a  significant  positive  correlation  be¬ 
tween  SeaWiFS  surface  daily  integrated  PAR  and  LUMCON  Lake 
Pontchartrain  surface  daily  integrated  PAR  (Fig.  2,  slope  =  0.91, 
R2  =  0.72,  p<0.01,  N  =  631).  Observed  and  predicted  values  of  sur¬ 
face  PAR  agreed  well  and  explained  90%  of  the  variance  of  the  valida¬ 
tion  data  set  (data  not  shown,  slope  =  1.10,  R2  =  0.90,  p<0.01, 
N  =  626).  Therefore  the  regression  coefficients  from  Fig.  2  were 
used  to  correct  satellite-derived  daily  integrated  PAR  for  LCS  analysis. 
The  regression  of  Gulf  Ecology  Division  KdPAR  versus  1/SD  (Fig.  3) 
resulted  in  a  slope  of  1.44  (R2  =  0.94,  p<0.01).  Of  the  1115  in  situ 
measurements  from  LUMCON  there  were  137  match-ups  with  Sea¬ 
WiFS  on  the  same  day.  Of  the  374  in  situ  measurements  from  Gulf 


(E  m'2  d*') 

Fig.  2.  Model  II  geometric  mean  linear  regression  of  daily  integrated  surface  PAR  from 
the  LUMCON  Lake  Pontchartrain  monitoring  station  and  SeaWiFS  flag  free  pixels  from 
2001  through  2007.  The  dashed  line  represents  the  one  to  one  fie 
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Ecology  Division  there  were  10  match-ups  with  SeaWiFS.  The  147 
paired  observations  from  LUMCON  and  Gulf  Ecology  Division  were 
randomly  divided  so  half  of  the  data  were  used  to  develop  regression 
coefficients  and  the  other  half  of  the  data  were  used  to  validate  the 
regression  analysis.  LUMCON  field  measurements  of  KdPAR  ranged 
from  0.05  to  14.4  m  1  with  a  mean  of  0.69  ±  1.01  S.D.  Gulf  Ecology 
Division  field  measurements  of  KdPAR  ranged  from  0.03  to  7.2  m_1 
with  a  mean  of  0.48  ±0.65  S.D.  The  regression  between  SeaWiFS 
Kd4qo  versus  combined  Gulf  Ecology  Division  KdPAR  and  LUMCON 
KdPAR  was  significant  (Fig.  4,  slope  =  1.07f  R2  =  0.65,  p<0.01, 
N  =  68).  Observed  and  predicted  values  of  diffuse  attenuation  coeffi¬ 
cients  were  in  reasonable  agreement  and  explained  68%  of  the  vari¬ 
ance  (data  not  shown,  slope  =1.07,  R2  =  0.68,  p<0.01,  N  =  79).  The 
regression  coefficients  from  Fig.  4  were  used  to  adjust  SeaWiFS 
Kd4qo  to  represent  KdPAR  hereafter  referred  to  as  Kd490  PAR. 

32.  Sparia/,  seasonal,  and  inrerannua/  patterns 

The  LCS  was  divided  into  five  discrete  regions  (Fig.  5),  to  assist  in 
spatial  descriptions.  The  percent  frequency  of  PAR  >  1%  at  the  SWI  for 
March,  April,  and  May,  from  1998  through  2007,  was  used  to  guide 
regional  divisions  because  this  was  the  time  of  year  when  hypoxia  be¬ 
gins  to  develop.  Regional  boundaries  were  drawn  based  on  the  visible 
spatial  patterns  in  Fig.  5.  Regions  were  ordered  in  decreasing  values  of 
average  Kd49o~PAR  and  increasing  average  water  depth.  The  outer  por¬ 
tion  of  the  domain  was  restricted  to  the  70  m  isobath.  Region  area, 
volume,  and  depth  were  reported  in  Table  1.  The  shallowest  mean 
depth  was  in  Region  1  (Rl)  and  the  euphoric  depth  reached  the 
SWI  <  5%  of  the  time  during  the  spring  months.  Region  2  (R2)  was 
just  south  of  Atchafalaya  Bay  and  the  euphoric  depth  reached  the 
SWI  >50%  of  the  time.  Region  3  (R3)  was  drawn  with  the  northern  ex¬ 
tent  along  the  10  m  isobath  and  the  euphoric  depth  reached  the  SWI 
between  0  and  40%  of  the  time  during  the  spring  season.  Region  4 
(R4)  was  located  adjacent  to  the  Mississippi  River,  Barataria  Bay, 
and  Terrebonne  Bay,  and  the  euphoric  depth  reached  the  SWI  <5% 
of  the  time.  Finally,  Region  5  (R5)  had  the  70  m  contour  along  its 
southern  most  extent.  Based  on  all  years  of  ship  observations,  R5 
was  characterized  by  the  highest  salinity  waters  and  the  euphoric 
depth  reached  the  SWI  >50%  of  the  time  during  March,  April,  and 
May. 

Kd4<x>~pAR  10-year  monthly  means  (Fig.  6A)  differed  significantly 
between  all  regions  for  every  month  of  the  year  (one-way  ANOVA, 
p<0.0001 ),  with  one  exception.  Kd4cxv  par  was  not  different  signifi¬ 
cantly  between  R3  and  R4  in  February  (one-way  ANOVA,  p  =  0.91). 
Kd49o-PAR  increased  significantly  from  May  to  June  (one-way 
ANOVA,  p< 0.001)  and  decreased  from  July  to  August  (p<  0.001)  in 
all  regions.  Rl  mean  Kd4qo  par  was  highest  followed  by  R2,  R3,  R4, 
and  R5  respectively.  Rl,  R2,  and  R3  means  Kd49a  PAR  peaked  in 


FlS- 3.  Model  11  geometric  mean  linear  regression  of  the  diffuse  attenuation  coefficient 
Kdp/m  determined  from  the  Gulf  Ecology  Division  CTD  PAR  data  and  inverse  of  Secchi 
depth. 


SeaWiFS  Kd4<K, 

<!»•') 

Fig.  4.  Model  11  geometric  mean  linear  regression  of  the  SeaWiFS  Kd^  diffuse  light  at¬ 
tenuation  coefficient  with  LUMCON  and  Gulf  Ecology  Division  Kd^.  The  dashed  line 
represents  the  one  to  one  fit. 

January  while  R4  and  R5  means  Kd4go_pAR  peaked  in  March  and  De¬ 
cember,  respectively.  Minimum  KcL^par  occurred  in  each  region 
during  the  months  of  August  or  September.  Seasonally,  KdUgo-PAR 
was  greater  in  the  winter  than  the  summer  in  all  regions.  Percent 
(Fig.  6B)  and  absolute  PAR  (Fig.  6C)  at  the  SWI  was  inversely  propor¬ 
tional  to  the  change  in  the  Kd49o~pAR-  Mean  percent  and  absolute  PAR 
at  the  SWI  was  lowest  in  the  fall  and  winter,  but  peaked  from  April 
through  September  with  a  decrease  in  July.  The  lowest  percent  and 
absolute  light  at  the  SWI  typically  occurred  in  December  and  January. 
SeaWiFS  daily  integrated  PAR  (Fig.  6B)  peaked  between  50  and 
60 Em  2d  1  during  June.  MARB  discharge  (Fig.  6D)  typically 
peaked  in  the  spring  and  was  at  a  minimum  from  August  through  Oc¬ 
tober.  Wind  speed  (Fig.  6D)  was  at  a  maximum  during  the  winter 
months  of  November  and  December  and  at  a  minimum  in  July  and 
August. 

MARB  discharge,  total  nitrogen  concentration,  and  wind  speed  in¬ 
fluences  on  Kd4<x)~PAR  were  used  to  assess  interannual  variability 
(Table  2)  and  shorter  time-scales  during  the  months  of  May,  June, 
and  July  (Table  3)  when  hypoxic  conditions  typically  proliferate. 
The  stepwise  multiple  linear  regression  analysis  on  the  interannual 
scale  (Table  2)  indicated  Kd490~PAR  was  correlated  with  wind  speed 
in  R5.  Both  wind  speed  and  discharge  were  significantly  related  to 
KCL490.  PAR  in  Rl,  R2,  and  R3.  Discharge  and  NOx  were  the  best  predic¬ 
tors  of  Kdjgo-PAR  in  R4.  Kd49o-pAR  during  May,  June,  and  July  (Table  3) 
was  mainly  related  to  Q  and  NOx.  For  these  months,  KdUgo^pAR  had  a 
significant  negative  correlation  with  wind  speed  in  R5  and  R4.  Corre¬ 
lations  among  independent  variables  (co-linearity)  may  produce  in¬ 
accurate  estimates  of  regression  coefficients,  variability,  and  p- 
values  (Kleinbaum  et  al.,  1986).  In  this  case,  the  lagged  Qand  lagged 
NOx  were  not  strongly  correlated,  but  literature  provided  minimal 
guidance  on  defining  what  level  of  correlation  constituted  a  problem. 
Thus,  only  one  Qand  one  NOx  value  was  used.  In  cases  where  multi¬ 
ple  lagged  variables  occurred,  only  the  variable  with  the  smallest  p- 
value  was  retained. 

The  SeaWiFS  derived  euphoric  depth  correlated  with  the  depth  at 
which  the  water  column  turned  hypoxic  at  <2mgL  1  (zHypoxic). 
explaining  51%  (p<0.001 )  of  the  variability  (Fig.  7A).  Shallow  eupho¬ 
ric  depths  occurred  where  there  was  shoaling  of  the  hypoxic  layer  to¬ 
ward  the  surface,  and  deeper  euphoric  depths  occurred  where  the 
hypoxic  layer  was  closer  to  the  SWI.  Mean  July  Kd49a  PAR  within  the 
maximal  extent  of  the  hypoxic  zone  (22,000  km2  reported  in  2002) 
was  correlated  to  the  normalized  size  of  the  July  hypoxia  area 
(Fig.  7B).  The  July  hypoxic  area  size  was  normalized  to  the  maximal 
extent  of  the  hypoxic  zone.  Normalized  hypoxic  area  increased  with 
an  increase  in  Kd49o-PAR-  Percent  surface  PAR  at  the  SWI  was  reported 
from  monthly  composites  between  1998  and  2007  (Fig.  8).  Mean  per¬ 
cent  PAR  increased  from  January  through  May  (Fig.  8A,  B,  C,  D,  and  E) 
within  all  regions.  R5  mean%  PAR  peaked  in  May  while  R2  and  R3 
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Fig.  5.  Northern  Gulf  of  Mexico  regions  based  the  spatial  pattern  of  the  percent  frequency  which  the  euphotic  layer  reaches  the  SW1  for  March.  April,  and  May  from  1998  through 
2007.  The  domain  of  interest  was  within  the  70  m  isobaths  The  70  and  10  m  isobaths  are  indicated  as  white  broken  lines. 


peaked  in  August  (Fig.  8H).  R1  and  R4  mean  %  PAR  peaked  in  Septem¬ 
ber  (Fig.  81).  Mean  percent  PAR  decreased  in  all  regions  from  October 
(Fig.  8J)  to  November  (Fig.  8K).  with  the  lowest  mean  percent  PAR 
reaching  the  SWI  in  December  (Fig.  8L).  Across  our  defined  domain 
46%  of  the  area  was  greater  than  0.1%  PAR  at  the  SWI  and  2%  of  the 
area  was  greater  than  1.0%  PAR  during  January  (Fig.  9).  Ninety-four 
percent  of  the  domain  was  greater  than  0.1%  PAR  and  66%  of  the 
area  was  greater  than  1.0%  PAR  during  May.  Finally,  during  July  90% 
of  the  domain  was  greater  than  0.1%  PAR  and  36%  of  the  area  was 
greater  than  1.0%  PAR. 


3.3.  Future  perspective 

Generally,  Kd^-PAR  decreased  under  nitrogen  reduction  and  in¬ 
creased  under  global  climate  change  scenarios  (Table  4).  R4  was  clos¬ 
est  in  proximity  to  the  Mississippi  River  and  had  the  largest  relative 
percent  change  during  all  scenarios.  R3  had  the  second  largest  per¬ 
cent  decrease  (  —  7%  and  — 12%)  under  the  30%  and  50%  nitrogen  re¬ 
duction  scenarios.  R2  and  R3  had  identical  increases  of  5  and  6% 
during  the  increased  nitrogen  and  discharge  climate  change  scenari¬ 
os.  R1  KcWpar  had  no  change  under  any  of  the  nitrogen  reduction 
scenarios  and  increased  6%  with  increased  discharge.  Finally,  R5 
Kd49o_PAR  had  no  change  under  the  increased  discharge  scenario,  in¬ 
creased  3%  with  increased  nitrogen,  and  decreased  5  and  8%  with 
30%  and  50%  nitrogen  reduction. 


Table  1 

Area,  volume,  and  depth  statistics  for  regions  1  through  5. 


Region 

Area 

(km2) 

Volume 
(xlO6  km3) 

Average  depth 

(m) 

Min  depth 

(m) 

Max  depth 

(m) 

R1 

1992 

9 

4.5 

1 

7 

R2 

3065 

19 

62 

1 

16 

R3 

9334 

135 

14.5 

3 

30 

R4 

4343 

123 

28.4 

7 

70 

R5 

35.084 

1293 

36.8 

7 

75 

Total 

53.818 

1579 

4.  Discussion 

4.1.  SeaWiFS  Kd490-PAR  and  PAR 

Differences  between  the  SeaWiFS  Kd49o  and  KdpAR  field  observa¬ 
tions  were  expected  because  KdPAR  was  the  integral  of  the  spectrum 
(400-700  nm)  and  Kd490  was  a  measure  at  a  single  waveband 
(Fig.  4).  The  PAR  spectrum  changed  with  depth  in  the  water  column 
due  to  different  absorption  rates  at  each  wavelength.  The  red  wave¬ 
lengths  attenuated  rapidly  with  depth  and  the  blue  wavelengths  pen¬ 
etrated  deeper  into  the  water  column  (Kirk,  1994).  Differences 
between  SeaWiFS  PAR  and  the  LUMCON  Lake  Pontchartrain  PAR  oc¬ 
curred  because  the  SeaWiFS  algorithm  was  calculated  from  a  single 
instantaneous  measurement  per  day  (Frouin  et  al.,  2003)  and  as¬ 
sumed  cloud  cover  did  not  change  throughout  the  day.  The  LUMCON 
PAR  was  influenced  by  changing  cloud  cover  throughout  the  day. 

Estimates  of  SeaWiFS  KdUgo  par  assumed  the  diffuse  light  attenua¬ 
tion  coefficient  was  constant  throughout  the  water  column.  Inspec¬ 
tion  of  Gulf  Ecology  Division  CTD  PAR  profiles  found  80%  of  vertical 
profiles  had  higher  diffuse  attenuation  in  the  surface  (unpublished 
data).  Lugo-Fernandez  et  al.  (2008)  also  reported  higher  diffuse  at¬ 
tenuation  in  surface  water  which  was  typically  over  clearer  Gulf  of 
Mexico  water.  Jahnke  etal.  (2008)  noted  ship-based  variations  in  dif¬ 
fuse  attenuation  measurements  were  biased  toward  calmer  sea  con¬ 
ditions  potentially  biasing  benthic  PAR  estimates.  The  first  optical 
depth  observed  by  SeaWiFS,  ranged  from  0.27  to  >  13  m.  The  shallow¬ 
est  first  optical  depth  occurred  in  R1  and  the  deepest  first  optical 
depth  occurred  in  R5.  The  SeaWiFS  data  record  provided  the  advan¬ 
tage  of  shelf-wide  spatial  coverage  and  long-term  temporal  measure¬ 
ments,  which  compensated  for  the  inability  of  the  satellite  to  measure 
Kd49o-PAR  throughout  the  water  column  to  the  SWI. 


42.  Sparia/.  seasonal ,  and  interannual  patterns  in  light  attenuation 

Kd49o_PAR  was  highest  in  the  late  fall,  winter,  and  early  spring  and 
lowest  in  the  summer.  This  trend  was  consistent  with  Gulf  Ecology 
Division  and  LUMCON  field  observations.  At  the  interannual  scale, 
patterns  in  Kd49o-PAR  co-varied  with  wind  speed  in  all  regions  except 
R4  (Table  2).  Previous  studies  of  diffuse  light  attenuation  on  the  LCS 
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Fig.  6.  Ten-year  mean  of  SeaWiFS  (A)  Kd49o  pa*  (m  ').  (B)  daily  integrated  surface 
PAR  {Em  2  d  \  black  dashed  line),  %  PAR  (solid  lines)  at  the  SW1,  and  (C)  absolute 
PAR  at  the  SW1  for  all  five  regions  each  month  of  the  year.  Error  bars  indicate 
2x  ±  SE.  Eight  year  mean  of  (D)  wind  speed  and  MARB  discharge. 


Table  2 

Multiple  regression  models  for  the  response  variable  SeaWiFS  Kd490  based  on  the  1998 
to  2007  monthly  averages  of  Kd,,*,  for  each  region  (Rl,  R2,  R3,  R4,  and  R5)  and  the  in¬ 
dependent  variables:  wind  was  the  monthly  average  wind  speed  (ms  ’)  for  each  re¬ 
gion,  Q  was  the  monthly  average  combined  discharge  (m3  s  1 )  of  the  Mississippi  and 
Atchafalaya  Rivers  and  Q,1  and  (12  were  the  one  and  two- month  lagged  average  dis¬ 
charges,  NOx  was  the  flow-weighted  mean  monthly  nitrate  concentration  for  the  com¬ 
bined  Mississippi  and  Atchafalaya  discharge  and  NOx.l  and  NOxJ  were  the  one  and 
two- month  lagged  average  nitrate  concentrations,  df  are  degrees  of  freedom.  p  are 
the  regression  coefficients.  Partial  p-values  are  bold- italicized  when  significant  at 
a  =  0.05. 


Region 

p-value 

R2 

df 

Variable 

Beta 

Partial  p-value 

Intercept 

Rl 

1.1E-16 

0.56 

90 

Wind 

Q 

0.09 

2.22E-05 

<00001 

<00001 

0.74 

R2 

2.2E-14 

05 

90 

Wind 

Q 

0.12 

5.05E-06 

<0.0001 

0.01 

-0.008 

R3 

c  IE-16 

0.62 

90 

Wind 

Q 

0.07 

2.77E-06 

<0.0001 

0.003 

0.01 

R4 

<  IE-16 

0.57 

90 

Q 

NOx.1 

1.12E-05 

5.76E-04 

<0.0007 

0.03 

0.24 

R5 

5.8E-13 

0.43 

94 

Wind 

0.026 

<00001 

0.02 

Table  3 

Multiple  regression  models  for  the  response  variable  SeaWiFS  Kd490  based  on  the  1998 
to  2007  May,  June,  and  July  averages  of  Kd^o  for  each  region  {Rl,  R2,  R3,  R4.  and  R5) 
and  the  independent  variables:  Wind  was  the  monthly  average  wind  speed  (ms"1) 
for  each  region.  Q  was  the  monthly  average  combined  discharge  (m3  s" ')  of  the  Mis¬ 
sissippi  and  Atchafalaya  Rivers  and  Q,1  and  Q,2  were  the  one  and  two-month  lagged 
average  discharges,  NOx  was  the  flow-weighted  mean  monthly  nitrate  concentration 
for  the  combined  Mississippi  and  Atchafalaya  discharge  and  NOx.1  and  NOx.2  were 
the  one  and  two-month  lagged  average  nitrate  concentrations,  df  are  degrees  of  free¬ 
dom.  p  are  the  regression  coefficients.  Partial  p-values  are  bold-italicized  when  signif¬ 
icant  at  a  =  0.05. 


Region 

p-value 

R2 

df 

Variable 

Beta 

Partial  p-value 

Intercept 

Rl 

0.015 

0.33 

21 

Wind 

-0.08 

0.22 

1.63 

Q 

2.18E-05 

0.005 

R2 

0.026 

024 

27 

Q 

828E-06 

0.07 

0.39 

NOx.2 

1.49E-03 

0.04 

R3 

2.0E-04 

0.61 

20 

Wind 

0.02 

024 

0.12 

Q 

6.12E-06 

00008 

NOx.2 

9.86E-04 

0.002 

R4 

3.0E-04 

0.6 

20 

Wind 

-0.09 

0.002 

0.5 

Q 

1.64E-05 

<00001 

NOx.2 

1.76E-03 

0.005 

R5 

0.009 

0.36 

21 

Wind 

-0.01 

0.02 

0.21 

NOx.2 

2.10E-04 

0.11 

were  limited  both  spatially  and  temporally  making  it  difficult  to  infer 
the  role  of  diffuse  light  attenuation  in  regulating  both  primary  pro¬ 
duction  and  hypoxia  across  the  shelf.  Diffuse  light  attenuation  coeffi¬ 
cients  reported  previously  (D’Sa  &  Miller,  2003;  Lehrter  et  al„  2009; 
Lugo- Fernandez  et  al„  2008)  were  similar  to  the  SeaWiFS  Kd490  par 
values  reported  in  this  study  throughout  regions  2,  3,  4,  and  5 
(Fig.  6).  Only  Region  1 ,  the  shallow  region  adjacent  to  the  Atchafalaya, 
had  higher  SeaWiFs  Kd49o_pAR  (2.64  m'1)  than  reported  previously. 

A  conceptual  explanation  of  the  possible  dynamics  was  developed, 
based  on  results  presented  here  and  previously  published  reports  that 


1998-2007 


Fig.  7.  Model  II  linear  regressions  of  (A)  euphotic  depth  estimated  from  SeaWiFS 
against  the  depth  at  which  the  water  column  turned  hypoxic  (zHypo*r.  <2.0  mg  L  '). 
and  (B)  mean  July  KcW  pa*  within  the  maximal  extent  of  the  hypoxic  area  against 
the  normalized  hypoxic  area  reported  for  the  month  of  July. 
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Fig.  8.  Ten-year  monthly  average  composites  ofSW]  percent  PAR  for  (A)  January.  (B)  February.  (C)  March.  (D)  April,  (E)  May.  (F)June,  (C)  July.  (H)  August.  (I)  September.  (J)  October, 
(K)  November,  and  (L)  December.  A  gray  land  mask  covers  coastal  waters,  such  as  the  Atchafalaya,  where  SeaWiFS  was  not  validated  with  field  observations. 


measured  the  individual  light  attenuating  components.  This  study  did 
not  quantify  individual  light  attenuating  components.  Diffuse  light  at¬ 
tenuation  is  the  sum  of  scattering  and  light  absorption  by 


Fig.  9.  Ten-year  monthly  averages  for  every  other  month  within  the  LCS  70  m  isobath 
that  has  less  than  or  equal  to  the  %  PAR  at  the  SW1. 


phytoplankton  pigments,  detritus,  and  colored  dissolved  organic  mat¬ 
ter  (CDOM).  Regions  in  close  proximity  to  the  MARB  discharge,  R1 
and  R4,  were  the  locations  where  Kc^-par  was  most  strongly  influ¬ 
enced  by  river  discharge  and  NOx  loading  (Tables  2  and  3).  Nitrogen 
loading  was  the  dominant  variable  controlling  phytoplankton  vari¬ 
ability  and  discharge  was  the  dominant  variable  controlling  detritus 
and  CDOM  variability  near  the  river  output  locations  (Green  & 
Gould,  2008;  Green  et  aL,  2008). 


Table  4 

Average  Kd  for  the  May  to  July  time  period  for  each  region.  The  30%  NOx  reduction,  50% 
NOx  reduction,  20%  NOx  increase,  and  20%  Q  increase  are  scenarios  that  were  evaluat¬ 
ed  as  to  their  effect  on  Kd.  Based  on  these  scenarios,  for  each  region  the  Kd  was  pre¬ 
dicted  using  the  multiple  regression  models.  Shown  are  the  average  percentage 
changes  in  predicted  Kd  from  the  average  Kd. 


Average 

Kd 

30%  NOx 
reduction 

50%  NOx 
reduction 

20%  NOx 
increase 

20%  Q 
Increase 

Region  1 

1.65 

0% 

0% 

0% 

6% 

Region  2 

0.72 

-6% 

-11% 

5% 

6% 

Region  3 

0.4S 

- 7 % 

-12% 

5% 

6% 

Region  4 

052 

-11% 

-18% 

7% 

13% 

Region  S 

0.1S 

-5% 

—  B% 

3% 

0% 
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Nutrient  and  discharge  rates  peaked  in  March,  April  and  May,  The 
increased  nutrient  availability  and  calm  winds  in  during  June,  July  and 
August  allowed  for  water  column  stabilization  and  phytoplankton  bio¬ 
mass  to  accumulate,  therefore  absorbing  more  PAR,  thus  increasing 
KdU9o-pAR  during  the  summer  months.  While  phytoplankton  absorption 
contributed  to  KcW  par  along  the  LCS  coastal  waters,  KdUqo-PAR  along 
the  coast  in  R1,  R2,  and  R4  was  dominated  by  CDOM  absorption, 
which  decreased  with  increasing  salinity  (D’Sa  &  Miller,  2003;  Green 
et  al„  2008).  The  dominant  CDOM  absorption  signature  and  influence 
of  river  discharge  in  MARB  proximal  regions  suggested  KcW  par  was 
controlled  not  only  by  phytoplankton  responses  to  nutrients,  but  also 
to  increased  CDOM  concentrations. 

Regions  2, 3,  and  5  Kd490-PAR  were  more  positively  correlated  with 
wind  speed  from  October  through  March.  Wind  speed  was  the  stron¬ 
gest  predictor  for  detritus,  while  discharge  and  wind  direction  were 
predictors  of  CDOM  (Green  et  al.,  2008).  Green  et  al.  (2008)  examined 
the  detritus  and  CDOM  absorption  components  in  relation  to  wind 
and  discharge,  but  it  was  not  discussed  how  these  terms  contributed 
to  total  diffuse  light  attenuation.  Higher  wind  speeds  increased  verti¬ 
cal  mixing  and  detritus  retention  in  the  water  column  increasing 
Kdm.PAR  (Salisbury  et  al..  2004).  As  river  waters  spread  across  the 
LCS  it  delivers  detritus  that  rapidly  sinks  below  the  pycnocline.  Sum¬ 
mertime  low  wind  speed  and  solar  heating  enhanced  the  buoyancy  of 
plume  waters  and  detritus  was  maintained  in  the  surface  layers 
(Walker  et  al..  1994). 

Winds  were  predominantly  from  the  east  except  during  June,  July, 
and  August  when  they  switched  to  southwest.  The  east  winds  pro¬ 
duced  downwelling  conditions  along  the  LCS  and  the  dominant  cur¬ 
rent  direction  was  to  the  west.  Southwest  winds  promoted 
upwelling  conditions  (Hetland  &  DiMarco,  2008).  The  change  to  off¬ 
shore  currents  and  upwelling  conditions  may  transport  MARB  surface 
water  onto  the  LCS,  and  may  also  increase  the  phytoplankton  biomass 
and  pigment  absorption  in  areas  such  as  R3  and  R5. 

The  Kd49o-PAR  of  R3  was  similar  to  R4  (Fig.  6A),  which  was  proxi¬ 
mal  to  a  MARB  river  output.  These  two  regions  were  typical  locations 
of  summer  hypoxia  on  the  LCS.  The  distinction  between  the  two  re¬ 
gions  was  the  percent  surface  PAR  which  reached  the  SWl  due  to  a 
difference  in  bathymetry,  R4  and  R3  average  depths  were  28.4  and 
14.5  m  (Table  1).  During  non-summer  months  LCS  current  direction 
was  to  the  west  which  adverted  Mississippi  River  discharge  over  R3 
and  downwelling  conditions  limited  Atchafalaya  discharge  from  pen¬ 
etrating  further  onto  the  LCS.  As  the  summer  winds  shifted  the 
Atchafalaya  plume  was  transported  eastward  toward  R3  and  R5. 
This  pattern  was  supported  by  Fig.  8  where  the  SWI  in  R3  received 
less  than  0.1%  surface  PAR  (purple  and  black)  October  (Fig.  8J) 
through  March  (Fig.  8C)  and  less  than  1%  surface  PAR  (green)  April 
(Fig.  8D)  through  September  (Fig.  81). 

Further  south  on  the  LCS,  R5  PAR  at  the  SWI  was  greater  than  1% 
surface  irradiance  (yellow  and  red)  and  extended  to  91.5°W  in  May 
(Fig.  8E)  and  June  (Fig.  8F),  but  was  limited  west  of  93.0°W  in  July 
(Fig.  8G)  and  never  regained  the  maximal  extent  by  the  end  of  the 
year  (Fig.  8H  through  8L).  Therefore  an  increase  in  Kd49o  par  during 
June  and  July  resulted  from  increased  phytoplankton  biomass  due  to 
maximal  nutrients,  discharge,  and  water  column  stabilization;  in  ad¬ 
dition  to  CDOM  and  detritus  advection  onto  the  LCS  due  to  changing 
winds  and  the  development  of  an  offshore  current. 

4.3.  Implications  for  primary  production  and  hypoxia 

Previous  studies  have  suggested  the  importance  of  light  in  regu¬ 
lating  sub-pycnocline  oxygen  dynamics  (Bierman  et  al.,  1994;  Lehrter 
et  al.,  2009).  This  study  supports  the  linkage  between  the  spatial  ex¬ 
tent  of  hypoxia,  both  vertically  (Fig.  7A)  and  horizontally  (Fig.  7B), 
and  were  correlated  to  Kd4qa.PAR  along  with  environmental  factors 
such  as  river  discharge,  nutrients,  and  wind  speed.  Changes  in 
Kd4qa  PAR  affect  sub-pycnocline  organic  matter  production  and  the 


primary  production  to  respiration  ratio  (P:R),  which  determined  the 
net  gain  or  loss  of  oxygen.  The  use  of  a  constant  threshold  of  1%  sur¬ 
face  light  for  phytoplankton  (Kirk,  1994)  or  0.1%  surface  light  for  MPB 
compensation  depths  (Cahoon,  1999)  was  a  suggested  indication  of 
light  levels  that  support  primary  production  and  hence,  oxygen  pro¬ 
duction.  Empirical  results  from  this  study  demonstrated  that  the  1% 
euphoric  depth  was  significantly  correlated  to  the  depth  of  hypoxia 
(Fig.  7).  Therefore,  a  euphoric  water  column  and  SWI  with  at  least 
0.1%  surface  irradiance  had  the  potential  for  net  oxygen  production 
(Gattuso  et  al.,  2006;  Glud  et  al.,  2002).  Thus,  low  KcU^o-par  may  pro¬ 
mote  sub-pycnocline  production  and  partially  offset  oxygen  deple¬ 
tion  by  respiration  below  the  pycnocline.  When  KcL^  PAR  increased, 
the  sub-pycnocline  P:R  ratio  decreased  and  the  net  drawdown  of  ox¬ 
ygen  proceeded  more  rapidly. 

The  increased  Kd49o-PAR  from  May  to  July  (Fig.  6)  coincided  with 
peak  stratification  strength,  thus  limiting  the  exchange  of  oxygen 
across  the  pycnocline.  There  were  two  implications  of  this  temporal 
pattern  as  regards  hypoxia.  First,  there  was  the  potential  for  large 
sub-pycnocline  primary  production  of  organic  matter  from  April  to 
May.  The  organic  matter  produced  in  these  waters  may  be  an  impor¬ 
tant  substrate  for  respiration  leading  to  the  development  of  hypoxia. 
Second,  the  increased  KcU^par  from  June  to  July  (Fig.  6)  would  result 
in  decreased  sub-pycnocline  primary  production.  Thus,  variation  in 
KdUgo  Par  (Figs.  7  and  8)  may  shift  the  sub-pycnocline  metabolism 
from  net  oxygen  production  to  net  oxygen  loss.  Eutrophication  and 
physical  processes  such  as  stratification  and  wind  mixing  events  pri¬ 
marily  control  the  extent  of  hypoxia,  while  light  may  partially  control 
the  P:R  ratio  in  the  sub-pycnocline  water  column  and  at  the  SWI.  In 
this  manner,  light  may  play  a  critical  role  in  structuring  the  ratio  of 
net  community  photosynthesis  and  respiration,  which  in  turn  influ¬ 
ences  the  timing  and  extent  of  hypoxic  water  mass  formation. 

4.4.  Future  perspective 

Nutrient  management  scenarios  of  30  and  50%  reductions  in  nitro¬ 
gen  had  the  greatest  impact  on  R4  and  R3.  It  could  be  expected  that 
R4  Kd490-PAR  would  have  the  greatest  percentage  change  due  to  the 
close  proximity  of  the  Mississippi  River,  but  it  was  surprising  that 
R3  was  the  second  largest  percentage  change.  R1  and  R2  were  closer 
in  proximity  to  the  Atchafalaya  River  discharge  source,  but  had  no  re¬ 
sponse  or  smaller  decreases  in  Kd4qo_PAR  relative  to  R3.  Global  climate 
change  scenarios  (Table  4)  reported  by  Miller  and  Russell  (1992)  and 
Wolock  and  McCabe  (1999)  suggested  a  20%  increase  in  MARB  dis¬ 
charge.  Increased  discharge  and  nitrogen  from  fertilizer  utilization 
also  increased  nitrogen  export  to  the  LCS  (Donner  &  Scavia,  2007; 
Donner  et  al.,  2002),  MARB  discharge  was  suggested  to  increase 
over  the  last  5  decades  (Dai  etal.,  2009),  and  forecasted  general  circu¬ 
lation  model  results  suggested  discharge  will  continue  to  increase 
(Nijssen  et  al.,  2001 ).  Increased  discharge  may  increase  diffuse  atten¬ 
uation  by  increasing  the  supply  of  nitrogen,  phytoplankton  biomass, 
CDOM,  and  detritus.  A  decrease  in  available  sub-pycnocline  light 
may  subsequently  change  the  balance  of  P:R  and  art  in  a  synergistic 
fashion  to  increase  the  area  and  severity  of  hypoxia  in  R2,  R3,  and  R4. 

5.  Conclusions 

The  SeaWiFS  daily  integrated  surface  PAR  was  validated  and  the 
KdUgo  diffuse  attenuation  product  was  adjusted  using  the  slope  and 
intercept  of  the  regression  against  field  observations  of  KdPAR.  At  the 
interannual  scale,  monthly  patterns  in  KcL^-par  were  strongly  corre¬ 
lated  with  wind  speed.  However,  Kd49a,PAR  was  related  to  MARB  dis¬ 
charge  and  nutrients  during  the  late  spring  and  summer.  Increases  in 
KdU9o„PAR  during  June  and  July  were  likely  the  result  of  increased  phy¬ 
toplankton  biomass  due  to  maximal  nutrients,  discharge,  and  water 
column  stabilization.  In  addition,  seasonal  dynamics  of  CDOM  and  de¬ 
tritus  may  vary  with  changes  in  wind  patterns.  There  was  potential 
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for  sub-pycnocline  primary  production  during  April  to  May.  Kd490_pAR 
increased  from  May  to  July,  which  may  shift  the  sub-pycnocline  me¬ 
tabolism  from  net  autotrophic  to  net  heterotrophic,  thus  potentially 
influencing  the  onset  of  hypoxia.  The  availability  of  light  is  one  of 
many  environmental  variables  that  interact  to  cause  hypoxia.  This 
paper  provides,  for  the  first  time,  an  analysis  of  light  availability 
across  the  large  area  of  the  LCS.  Our  study  was  also  the  first  to  show 
that,  in  addition  other  factors  such  as  stratification,  river  discharge, 
and  nutrients,  the  vertical  distribution  of  hypoxia  in  the  water  col¬ 
umn  was  related  to  the  euphotic  depth  and  the  horizontal  extent  of 
hypoxia  was  also  related  to  the  change  in  Kd490_PAR.  Finally,  global  cli¬ 
mate  change  scenarios  suggest  Kd490,pAR  may  increase  across  the  LCS 
potentially  increasing  the  hypoxic  area  and  shoaling  the  depth  at 
which  the  water  column  turned  hypoxic. 

Acknowledgments 

This  study  has  been  funded  in  part  by  the  U.  S.  Environmental  Pro¬ 
tection  Agency.  It  has  been  reviewed  and  approved  for  publication  by 
the  National  Health  and  Environmental  Effects  Research  Laboratory 
but  contents  are  solely  views  of  the  authors.  Mention  of  trade 
names  or  commercial  products  does  not  constitute  endorsement  by 
the  US  EPA.  This  is  contribution  number  1375  from  the  Gulf  Ecology 
Division.  George  Smith  provided  G1S  support.  We  thank  Don  Redalje 
at  the  University  of  Southern  Mississippi  Department  of  Marine  Sci¬ 
ence  for  his  comments  and  the  anonymous  reviewers  for  their 
suggestions. 

References 

Amone,  R.  A.,  Martinolich.  P.  M.,  Could.  R.  W„  5ydor,  M„  5tumpf,  R.  P„  A  Ladner.  5. 

(1998).  Coastal  optical  properties  using  SeaWiFS.  Oceon  optics  XIV.  HA:  Kona. 
Bailey.  5.  W.,  A  Werdell.  P.  J.  (2006).  A  multi-sensor  approach  for  the  orbit  validation  of 
ocean  color  satellite  data  products.  Remote  Sensing  of  Environment.  102,  12-23. 
Baith,  1C  Lindsay.  R.,  Fu,  G,,  &  McClain.  C  R.  (2001).  SeaDAS:  Data  analysis  system  de¬ 
veloped  for  ocean  color  satellite  sensors.  EOS.  Transoctions,  American  Geophysicol 
Onion,  82. 

8ierman,  V.  J..  Hinz,  5.  C.,  Wiseman.  W.  J.,  Rabalais,  N.  N.,  &  Turner,  R.  E.  ( 1994).  A  pre¬ 
liminary  mass  balance  model  of  primary  productivity  and  dissolved  oxygen  in  the 
Mississippi  River  plume/inner  Gulf  5helf  region.  Estuories,  17,  886-899. 

Cahoon,  L  B.  ( 1999).  The  role  of  benthic  microalgae  neritic  ecosystems.  Oceonography 
ond  Morine  Biology,  37.  47-86. 

Dai,  A.  G.,  Qian,  T.  T.,  Trenberth,  K.  L,  &  Milliman,  J.  D.  (2009).  Changes  in  continental 
freshwater  discharge  from  1948  to  2004.  Joumol  of  Climote,  22.  2773-2792. 

Dale.  V.  H..  Kling,  C.  L.  Meyer.  J.  L.  Sanders.  J..  Stallworth,  H.,  Armitage,  T,  et  al.  (2010). 

Hypoxio  in  the  Northern  Gu//a/ Mexico.  New  York:  Springer. 

Donner,  5.  D.,  Coe,  M.  T„  Lenters.  J.  D..  Twine,  T.  L,  A  Foley.  J.  A.  (2002).  Modeling  the 
impact  of  hydrological  changes  on  nitrate  transport  in  the  Mississippi  River  Basin 
from  1955  to  1994.  Global  Biogeochemical  Cycles.  16.  doi:  1 0.1 029/200 1GBOO 1396. 
Donner,  5.  D„  &  Scavia,  D.  (2007).  How  climate  controls  the  flax  of  nitrogen  by  the  Mis¬ 
sissippi  River  and  the  development  of  hypoxia  in  the  Gulf  of  Mexico.  Limnology  ond 
Oceonography  fUmnoL  Oceanogr.J,  52.  856-861. 

Dortch.  Q„  Rabalais,  N.  N„  Turner,  R.  EL,  &  Rowe.  G.  T.  (1994).  Respiration  rates  and  hyp¬ 
oxia  on  the  Louisiana  shelf.  Estuories.  1 7,  862-872. 

D*Sa,  E.  J.,  A  Miller.  R.  L  (2003).  8io-optical  properties  in  waters  influenced  by  the  Missis¬ 
sippi  River  during  low  flow  conditions.  Remote  Sensirtg  of  Environment,  84, 538-549. 
Eadie,  B.J„  McKee.  B.  A.  Lansing,  M.  8.,  Robbins,;  A,  Metz,  5„  ATrefry.J.  H.  (1994).  Re¬ 
cords  of  nutrient-enhanced  coastal  productivity  in  sediments  from  the  Louisiana 
continental  shelf.  Estuones,  17, 754-765. 

Feldman.  G.  C,  &  McClain.  C  R.  (2008).  Ocean  Color  Web,  SeaWiFS  reprocessing  5.2.  In 
N.  Kuring,  &  5.  W.  Bailey  (Eds.),  NASA  Goddord  Spoce  Flight  Center  http://oceancolor. 
gsfc.nasa.gov 

Frouin,  R^  Franz,  8„  &  Werdall,  P.  J.  (2003).  The  SeaWiFS  PAR  product  In  S.  8.  Hooker,  & 
E.  R.  Firestone  (Eds.),  SeoWiES  postiounch  technicol  report  series  (pp.  46-50).  (Algo¬ 
rithm  Updates  for  the  Fourth  SeaWiFS  Data  Reprocessing:  NASA  Technical 
Memorandum). 

Frouin.  R.,  Franz,  8.  A..  &  Werdell,  P.  J.  (2003).  The  SeaWiFS  PAR  product  In  S.  8.  Hooker, 
&  E.  R.  Firestone  (Eds.).  SeaWiFS  postiounch  technicol  report  series  (pp.  46-50).  (Al¬ 
gorithm  Updates  for  the  Fourth  SeaWiFS  Data  Reprocessing). 

Gattuso,  J.  P.,  Gentili.  B..  Duarte.  C  M..  Kleypas,  J.  A.  Middelburg,  J.  J„  &  Antoine,  D. 
(2006).  Light  availability  in  the  coastal  ocean:  Impact  on  the  distribution  of  benthic 
photosynthetic  organisms  and  their  contribution  to  primary  production.  Biogeos¬ 
ciences,  3,  489-513. 

Glud,  R.  N..  Kuhl,  M..  Wenzhofer.  F„  &  Rysgaard,  S.  (2002).  Benthic  diatoms  of  a  high 
Arctic  Qord  (Young  Sound.  NE  Greenland):  Importance  for  ecosystem  primary  pro¬ 
duction.  Morine  Ecology-Progress  Series.  23B,  15-29. 


Goolsby,  D.  A,  Battaglin,  W.  A,  Aulenbach,  B.  T„  &  Hooper.  R.  P.  (2001 ).  Nitrogen  input 
to  the  Gulf  of  Mexico,  /oumol  of  Environmental  Quolify,  30,  329-336. 

Green,  R.  E^  A  Gould.  R.  W,  (2008).  A  predictive  model  for  satellite-derived  phyto¬ 
plankton  absorption  over  the  Louisiana  shelf  hypoxic  zone:  Effects  of  nutrients 
and  physical  forcing,  Joumol  of  Geophysicol  Reseorrh,  113.  dol: 
10.1029  2007JC004 594. 

Green.  R.  E^  Gould,  R.  W„  &  Ko,  D.  5.  (2008).  Statistical  models  for  sediment/detritus 
and  dissolved  absorption  in  coastal  waters  of  the  northern  Gulf  of  Mexico.  Conti¬ 
nental  5hel/ Research.  28, 1273-1285. 

Greene.  R.  E,  Lehrter.  J.  C,  &  Hagy,  J.  D.  (2009).  Multiple  regression  models  for  hind  ■ 
casting  and  forecasting  midsummer  hypoxia  in  the  Gulf  of  Mexico.  Ecoiogicol  Appli- 
rations,  19.  1161-1175. 

Hetland,  R.  D.,  &  DiMarco.  5.  F.  (2008).  How  does  the  character  of  oxygen  demand  con¬ 
trol  the  structure  of  hypoxia  on  the  Texas- Louisiana  continental  shelf?  Joumol  of 
Marine  5ys terns,  70, 49-62. 

Jahnke,  R.  A,  Nelson.  J.  R^  Richards.  M.  E..  Robertson,  C  Y„  Rao,  A  M.  F.,  &  Jahnke,  D.  8. 
(2008).  Benthic  primary  productivity  on  the  Georgia  midcontinental  shelf:  Benthic 
flux  measurements  and  high-resolution,  continuous  in  situ  PAR  records.  Joumol  of 
Geophysicol  Research,  1 13.  doi:  10.1029  2008JC004745. 

JPL  (2001).  QuickScat  science  data  product  user's  manual:  Overview  and  geophysical 
data  products.  PubL  D-l  8053.  Pasadena,  Colifowio 

Justic,  D„  Rabalais,  N.  N.,  &  Turner,  R.  E.  (1996).  Effects  of  climate  change  on  hypoxia  in 
coastal  waters:  Doubled  CO2  scenario  for  the  northern  Gulf  of  Mexico.  Limnology 
ond  Oceanography,  41,  992-1003. 

Kirk,  J.  T.  0.  (1994).  Light  ond  photosynthesis  in  oquoric  ecosystems.  New  York:  Cam¬ 
bridge  University  Press. 

Kleinbaum,  D.  G.,  Kupper,  L  L,  A  Muller,  K.  E.  (1986;  Applied  regression  onolysis  ond 
other  multivariable  methods.  Boston:  PW5-Kent  Publishing, 

Laws,  E.  A,  &  Archie.J.  W.  (1981 ).  Appropriate  use  of  regression  analysis  in  marine  bi¬ 
ology.  Morine  Biology.  65. 13-16. 

Lee.  Z.  -P.,  Darecki.  M..  Carder,  K.  L,  Davis,  C  CL,  5tram$ki,  D.,  &  Rhea,  W.  J.  (2005).  Dif¬ 
fuse  attenuation  coefficient  of  downwelling  irradiance:  An  evaluation  of  remote 
sensing  methods.  Joumol  of  GeophysiccA  Rescorch,  1 10,  doi:  10. 1029/2 004JC0025 73, 

Lehrter.  J.  C  Kurtz,;  C,  A  Murrell.  M.  C  (2009).  Interactions  between  Mississippi  River 
inputs,  light,  phytoplankton  biomass,  and  phytoplankton  production  on  the  Louisi¬ 
ana  Continental  Shelf.  Continental  Shelf  Research,  29. 1861-1872. 

Lohrenz,  5.  E,  Fahnenstiel,  G.  L,  Kirkpatrick,  G.  J„  Carroll,  C.  L,  A  Kelly,  K.  A  (1999).  Mi¬ 
crophotometric  assessment  of  spectral  absorption  and  its  potential  application  for 
characterization  of  harmful  algal  species.  Joumol  0/ Phycalogy,  35.  1438-1446. 

Lohrenz,  5.  E,  Fahnenstiel,  G.  A,  Redalje.  D.  G.,  Lang.  G.  A^  Chen.  X.,  A  Dagg,  M.  J.  ( 1 997). 
Variations  in  primary  production  of  northern  Gulf  of  Mexico  continental  shelf  wa¬ 
ters  linked  to  nutrient  inputs  from  the  Mississippi  River.  Morine  Ecology  Progress 
Series,  155,  45-54. 

Lohrenz,  5.  E,  Redalje,  D.  G„  Cai,  W.  -J.,  Acker,  J.,  A  Dagg,  M.  J.  (2008).  A  retrospective 
analysis  of  nutrients  and  phytoplankton  productivity  In  the  Mississippi  River 
plume.  Continental  Shel/Researrh,  28,  1466-1475. 

Lugo-Femandez,  A,  Gravois,  M.,  A  Montgomery,  T.  (2008).  Analysis  of  secchi  depths 
and  light  attenuation  coefficients  in  the  Louisiana-Texas  shelf,  northern  Gulf  of 
Mexico.  Gulf  0/ Mexico  Science,  1, 14-27. 

Martinolich,  P.  M.  (2005).  Automated  processing  system  user’s  guide  v.3.4.,  Naval  Re¬ 
search  Laboratory.  Washington.  D.C  www7333,nr!sse.na vy.mil/docs  aps  v3.4/ 
user/aps 

Miller,;  R.,  A  Russell,  G.  L  (1992).  The  impact  of  global  warming  on  river  discharge. 
Joumol  0/ Geophysical  Reseorch.  97.  2757-2764. 

Mississippi  River/Gulf  of  Mexico  Watershed  Nutrient  Task  Force  (2001 ).  Action  plan  for 
reducing,  mitigating,  and  controlling  hypoxia  in  the  Northern  Gulf  of  Mexico. 
Washington,  DC 

Mississippi  River/Gulf  of  Mexico  Watershed  Nutrient  Task  Force  (2008).  Gulf  hypoxia 
action  plan  2008  for  reducing,  mitigating,  and  controlling  hypoxia  in  the  Northern 
Gulf  of  Mexico  and  improving  water  quality  in  the  Mississippi  River  Basin.  Wash¬ 
ington,  DC 

Nijssen,  B.,  O’Donnell,  G.  M.,  Hamlet,  AF,A  Lettenmaiet,  D.  P.  (2001).  Hydrologic  sen¬ 
sitivity  of  global  rivers  to  climate  change.  Climote  Chonge,  50,  143-175. 

Rabalais.  N.  N.,  Atilla,  N„  Normandeau.  C„  A  Turner,  R.  E  (2004).  Ecosystem  history  of 
Mississippi  River-Influenced  continental  shelf  revealed  through  preserved  phyto¬ 
plankton  pigments.  Morine  Pollution  Bulletin,  49,  537-547. 

Rabalais,  N.  N„  Turner.  R.  E,  Gupta,  B.  K.  5.,  Platon,  E^  A  Parsons,  M.  L  (2007).  Sediments 
tell  the  history  of  eutrophication  and  hypoxia  in  the  northern  Gulf  of  Mexico.  Eco¬ 
logical  Applications,  17.  SI  29-S 1 43. 

Redalje,  D.  G.,  Lohrenz,  5.  E,  A  Fahnenstiel,  G.  L  (1994).  The  relationship  between  pri¬ 
mary  production  and  the  vertical  export  of  particulate  organic  matter  In  a  river  im¬ 
pacted  coastal  ecosystem.  Estuories,  17,  829-838. 

Salisbury.].  E,  Campbell,;  W„  Linder,  E,  Meeker,  L  D„  Muller-Karger,  F.E,  A  Vorosmarty, 
CJ.  (2004).  On  the  seasonal  correlation  of  surface  particle  fields  with  wind  stress  and 
Mississippi  discharge  in  the  northern  Gulf  of  Mexico.  Deep-Seo  Research  U,  51, 
1187-1203. 

Scavia,  D„  Rabalais,  N.  N.,  Turner,  R,  E,  Justic,  D.,  A  Wiseman,  W.  J„  Jr.  (2003).  Predicting 
the  response  of  Gulf  of  Mexico  hypoxia  to  variations  in  Mississippi  River  nitrogen 
load.  Limnology  ond  Oceonography,  4B,  951-956. 

5tumpf,  R.  P„  Amone,  R.  A,  Gould,  R.  W„  Martinolich,  P.,  A  RansibrahmanakuL  V.  ( 2003).  A 
partially-coupled  ocean-atmosphere  model  for  retrieval  of  water-leaving  radiance 
from  SeaWiFS  in  coastal  waters.  In  5. 8.  Hooker  (Ed.),  A(gorir/im  Updates  /or  the  fourth 
SeaWiFS  Data  Reprocessing  {pp.  51-59). :  NASA  Technical  Memo. 

Turner,  R.  E,  Milan,  C  S.,  A  Rabalais,  N.  N.  (2004).  A  retrospective  analysis  of  trace 
metals,  C,  N,  and  diatom  remnants  in  sediments  from  the  Mississippi  River  delta 
shelf.  Morine  Pollution  8ulletin,  49, 548-556. 


BA  Schoeffer  et  oL  /  Remote  Sensing  of  Environment  1 15  (201 1 )  3745-3757 


3757 


Turner,  R.  L,  &  Rahalais,  N.  N.  ( 1991).  Changes  in  the  Mississippi  River  quality  in  this 
century.  BioSd,  41, 140-147. 

Walker,  N.  D„  Fargion,  C.  S„  Rouse,  L  J„  &  Biggs,  D.  C  ( 1 994).  The  great  flood  of  summer 
1993:  Mississippi  River  discharge  studied.  EOS.  Transactions,  75(409),  414-415. 

Wolock.  D.  M„  &  McCabe,  C.  J.  (1999).  Estimates  of  runoff  using  water-balance  and  at¬ 
mospheric  general  circulation  models.  Journo!  of  the  American  Woter  Resources  As¬ 
sociation,  35, 1341-1350. 


Wysocki,  L  A,  Bianchi,  T.  S„  Powell,  R.T..&  Reuss,  N.  (2006).  Spatial  variability  in  the 
coupling  of  organic  carbon,  nutrients,  and  phytoplankton  pigments  in  surface  wa¬ 
ters  and  sediments  of  the  Mississippi  River  plume.  Esruorfne,  Coostoi  ond  She!/  Sci¬ 
ence,  69, 47-63. 


